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Observations of spiral galaxies strongly support a one-to-one analytical relation between the 
inferred gravity of dark matter at any radius and the enclosed baryonic mass. It is baffling that 
baryons manage to settle the dark matter gravitational potential in such a precise way, leaving no 
"messy" fingerprints of the merging events and "gastrophysical" feedbacks expected in the history 
of a galaxy in a concordance Universe. This correlation of gravity with baryonic mass can be 
interpreted from several non-standard angles, especially as a modification of gravity called TeVeS, 
in which no galactic dark matter is needed. In this theory, the baryon-gravity relation is captured by 
the dieletric-like function /i of Modified Newtonian Dynamics (MOND), controlling the transition 
from 1/r^ attraction in the strong gravity regime to 1/r attraction in the weak regime. Here, we 
study this /x-function in detail. We investigate the observational constraints upon it from fitting 
galaxy rotation curves, unveiling the degeneracy between the stellar mass-to-light ratio and the 
/i-function as well as the importance of the sharpness of transition from the strong to weak gravity 
regimes. We also numerically address the efi'ects of non-spherical baryon geometry in the framework 
of non-linear TeVeS, and exhaustively examine how the function connects with the free function 
of that theory. In that regard, we exhibit the subtle effects and wide implications of renormalizing 
the gravitational constant. We finally present a discontinuity-free transition between quasi-static 
galaxies and the evolving Universe for the free function of TeVeS, inevitably leading to a return to 
attraction at very low accelerations in isolated galaxies. 

PACS numbers: 98.10.-(-z, 98.62.Dm, 95.35.-fd, 95.30.Sf 



I. INTRODUCTION 

Data on large scale structures point towards a Universe 
dominated by dark matter (DM) and dark energy [l|, Q ■ 
Nowadays, the dominant paradigm is that DM is actually 
made of non-baryonic weakly-interacting massive parti- 
cles, the so-called cold dark matter (CDM). However, on 
galaxy scales, the observations appear to be at variance 
with a sizeable list of GDM predictions. Measurements of 
non-circular motions in the Milky Way have shown that 
there is actually very little room for an axisymmetric 
distribution of dark matter inside the solar radius [l, |3| , 
where CDM simulations predict a cuspy density profile 
0. Even though CDM halos themselves could well be 
triaxial the non-circular motion of gas in the inner 
Milky Way is proven to be linked with non-axisymmetric 
baryonic features such as the bar and the spiral arms Q , 
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thus leaving no room for a cuspy triaxial CDM halo to 
explain them. External galaxies have also been used to 
compare the predicted cuspy CDM density profiles with 
the observations, in particular rotation curves of dwarf 
and spiral galaxies. Despite the presence of observational 
systematic effects [T], most of which are now under con- 
trol, observations point towards dark matter halos with 
a central constant density core [1, [13, [HI, [HI at odds 
with the CDM predictions. Another interesting problem 
faced by CDM on galactic scales is the overabundance 
of predicted satellite galaxies compared to the observed 
number in Milky Way-sized galaxies [13] ■ 

What is more, it is now well-documented that rotation 
curves suggest a correlation between the mass profiles 
of the baryonic matter (stars -I- gas) and dark matter 
[H) [I3| • Some rotation curves even display obvious 
features (bumps or wiggles) that are also clearly visible 
in the stellar or gas distribution [17|. A solution to all 
these problems, and especially the baryon-DM relation, 
could be some kind of new specific interaction between 
baryons and some exotic DM [l^. On the other hand. 
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it could indicate that, on galaxy scales, the observed dis- 
crepancy rather reflects a breakdown of Newtonian dy- 
namics in the ultra-weak field regime: this original pro- 
posal, dubbed MOND (Modified Newtonian Dynamics), 
was made more than twenty years ago in a series ofhall- 
mark papers [H, [l9|, 1131 • But as noted in, e.g., a 
modified gravity theory is in many ways equivalent to an 
exact relation between baryons and dark matter. When 
analyzing rotation curves, the two interpretations are 
nearly degenerate in the sense that the fraction of gravity 
unaccounted by the Newtonian gravity of the baryons at 
radius R is 

^ _ ^GNMb{R) ^ { ^"^^"^-^^"^^ in DM picture 

K^g |[l-^(5/ao)] in MOND picture 

(1) 

where Gat is the usual Newtonian gravitational constant, 
Mb{R) is the baryonic mass interior to the galactocentric 
radius i?, ^'^^^^^(R) (denotes the Newtonian gravity of the 
baryons (with a shape factor ^ ~ 1 to correct for non- 
spherical baryon distribution), and g = v^/Ris the total 
centripetal gravitational acceleration from measuring the 
circular speed v at radius R. The first equality introduces 
a spherical halo of dark matter with mass Mdm(-R) en- 
closed inside radius R (a flattening factor can also be 
used for the DM halo). The second equality (Milgrom's 
law) invokes a modifled dynamics with a dielectric-like 
[23 , [23I factor n < 1 to boost Newtonian gravity below a 
dividing acceleration scale qq whose value turns out to be 
of the same order as cHq . This fi factor is also called the 
MOND interpolating function /i(a;) (where x = g/ao)- 

The key difference between the usual dark matter and 
the MOND interpretation is the fact that MOND requires 
all galaxies to share the same /^-function for the same 
value of oq. If it is not the case, then the present-day 
success of MOND just tells us something phenomenolog- 
ical about DM. On the other hand, if it is indeed exactly 
the case, then either a to-be-discovered exotic relation 
does actually exist between DM (not the usual CDM) 
and baryons, or the DM paradigm is incorrect on galaxy 
scales. In the latter case, either gravity should be modi- 
fied [H, [2^ , which could lead to some subtle differences 
regarding e.g. dvnamical fricti on [26l I27I o r tides [28ll29j. 
or inertia should be modified [2ll ISOl. IsH. which would 
lead to very significant differences with the DM interpre- 
tation, especially in non-stationary situations 21]. 

The basis of the MOND paradigm J,8] is that /i is a 
function which runs smoothly from ^{x) = 1 at a; 3> 1 to 
/i(x) = a; at a; <C 1. Even though this leaves some free- 
dom for the exact shape of /i, it is already extremely sur- 
prising, from the DM point of view, that such a prescrip- 
tion did predict many aspects of the dynamics of galaxies 
such as the appearance of large discrepancies with New- 
tonian dynamics in low surface brightness galaxies, as 
well as the absence of discrepancies in the inner parts 
of high surface brightness ones. With this definition of 
the //-function, the fraction of gravity unaccounted by 
Newtonian gravity tends to zero in the high-acceleration 



limit, and one recovers in the low- acceleration limit the 
TuUy-Fisher law [H, [s^ , which tells us that the circu- 
lar velocity at the last measured radius in spiral galaxies 
must behave like M^^^^. Since the centripetal acceler- 
ation at that radius is typically lower than ~ 0.4ao in 
high surface brightness galaxies (and even lower in low 
surface brightness galaxies), the MOND circular veloc- 
ity at that radius is within at most ten percent of the 
asymptotic one, depending on the exact shape of the ^l- 
function, and on the stellar mass-to-light ratio. Indeed, 
the exact shape of the /i-function is still somewhat am- 
biguous in MOND, because there is a degeneracy between 
the choice of /j, and the stellar mass-to-light ratio when 
rotation curves are fitted with Eq. (1). While the exact 
asymptotic velocity will actually be different depending 
on the mass-to-light ratio, the shape of the /Lt-function will 
at the same time determine at what radius the rotation 
curve really becomes flat, and this can happen somewhat 
away from the last observed radius, thus explaining this 
degeneracy. 

Even though the most important aspect of MOND, 
confronting CDM with the most severe observational 
difficulties, remains its asymptotic behaviour described 
above, the study of the transition between the Newtonian 
regime and the asymptotic MOND regime, characterized 
by the shape of the /i-function, may also be important 
in the quest for a possible fundamental theory underpin- 
ning the MOND prescription. In this contribution, we 
thus investigate the observational constraints upon the 
/t-function in that transition zone, as well as the theoret- 
ical constraints upon it in the framework of Bekenstein's 
relativistic multifield theory of gravity [25j . 

II. OBSERVATIONAL CONSTRAINTS ON THE 
/t-FUNCTION 

A. The a-family and the standard function 

In this section, we examine the observational con- 
straints upon the /i-function from fitting galaxy rotation 
curves with Milgrom's law (the second equality of Eq. 1). 
As a matter of fact, we shall prove in Sect. Ill F that 
this empirical law yields an excellent approximation to 
Bekenstein's relativistic multifield theory of MOND [2^ 
in fitting rotation curves of spiral galaxies. 

Here, we explicitly show the degeneracy between the 
stellar mass-to-light ratio (-f the distance of the galaxy) 
and the shape of the /i-function, a degeneracy that pre- 
vents from fixing the /i-function in present-day MOND. 
Among a parametric family of //-functions, we are going 
to show that the best fits are obtained for the "simple" 
/i-function of [1]: 



and that, notwithstanding the M^,/L-fj, degeneracy 
(where M^/L is the stellar mass-to-light ratio), /i- 
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x = g/a(, 

FIG. 1: The /i-functions of Eq. (H) for a = 0, a = 0.5, a = 0.9, 
and a = 1, illsutrating that the transition from the deep- 
MOND [fj, ~ x) to the Newtonian (/x ~ 1) regime is sharper 
with increasing ct. 



TABLE I: Best-fit distances (d in Mpc) and stellar M/L ra- 
tios (in solar units) of the MOND rotation curve fits with the 
standard /i-function (Eq. |3} and the simple /i- function (Eq. [2] 
a = 1). Distances were left as free parameters within 30% of 
the distances adopted in the original papers; for NGC 7331 
the bulge M/L was considered as an additional free param- 
eter, with a best-fit value of 1.5 for both /^-functions. The 
degeneracy between those free parameters and the choice of 
the /i-function appears clearly (A'h/L playing the most im- 
portant role in this regard). 
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functions with too slow a transition from the Newto- 
nian to the deep-MOND regime are excluded by rotation 
curves data, while the "standard" /i- function [3j| , 



pi 



(3) 



with a faster transition than Eq. 12]), yield fits of compa- 
rable quality. 

The set of /z-functions we propose to test against ob- 
servations is called the a-family ^35i] : 



2x 



1 + (2 - a)x + ^(1 - axY + Ax' 



< a < 1. 



(4) 

In the a — \ case, one directly recovers the simple func- 
tion of Eq. In the a — Q case, it is a simple exercise 
to show that one recovers the /i-function of Bekenstein's 
toy-model [H, eq.(62)] 



^i{x) 



-1 + Vl + 4a; 



1 + 



(5) 



Fig. 1 displays those as a function of x. This is of 
course far from testing all the other possibilities, notably 
the free function proposed in [s^ , but it will give a good 
indication of the range of allowed /z-functions to be com- 
patible with observed rotation curves. 



B. Fitting galaxy rotation curves 

In order to test observationally the parametric a- 
family of /x-functions, we selected a sample made of two 



subsamples of rotation curves: the first subsample is 
made of nine rotation curves presented in (37| . This 
sample has been used to determine the value of the accel- 
eration constant ao — 1.21 x 10~® cm s~^ that is gener- 
ally used in MOND studies with the standard /i-function. 
The second subsample consists of the five galaxies pre- 
sented in @ . This way, we have a final sample of 14 galax- 
ies with high-quality, extended rotation curves. These 
galaxies span a large range of Hubble types and maximal 
velocities, from small irregular dwarf galaxies to large 
early-type spirals. 

The values of and a are unknown a priori, but have 
to be the same for all the galaxies. Therefore, we per- 
formed a global fit for these two parameters to the above- 
mentioned 14 rotation curves: a was allowed to vary 
from to 1, while the stellar M/L ratio and the distance 
(±30% the distances adopted in [sj and [§]) were left as 
parameters free to vary individually for each galaxy (see 
Table I). 

The best fits are obtained for 



a = 1 



+0.0 
-0.5' 



and 



ao 



l-35±o:42 X 10"^cms"^ 



(6) 



(7) 



where the errors are the one-sigma uncertainties from the 
statistics of the global fit; this value of ao is compati- 
ble with the result of [sj- For the unfavoured functions 
with a < 0.5, the transition from the deep-MOND to the 
Newtonian regime is slower than for a > 0.5 as shown on 
Fig. 1. 
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The results of the rotation curve fits using the best fit 
parameters a = 1 and ao = 1.35 x 10~^ cm s~^ are shown 
in Fig. [21 When comparing the fits of Fig. [2] to those 
performed with the standard MOND /i-function (EqO 
Fig. [3]), we see that the simple /i- function clearly fits the 
rotation curves as well as the standard one: the values 
of the fits are very similar. On average, we have: 



xLi/xL„da.d = 0.98 ±0.24. 



(8) 



Note however that three galaxies (NGC 2403, NGC 3198, 
and NGC 7331) also have Cepheid-based distances [ssl ]. 
While for NGC 2403 and NGC 7331 the best-fit distances 
(of both /x-functions) are within 10% of the Cepheid dis- 
tances, the best-fit distance of NGC 3198 is about a fac- 
tor 2 smaller than the Cepheid-based one. Imposing this 
value of the distance slightly decreases the best-fit value 
of ao but it doesn't change the best-fit value of a, nor 
does it change significantly the result of Eq. ([5]) , because 
the fits with both /i functions get worse for NGC 3198. 

The relatively large error bars on a are due to the de- 
generacy between the shape of /i and the stellar AI/L 
ratio. Still, it is striking that, notwithstanding this de- 
generacy, a = 1, corresponding to the simple relation 
of Eq. (21), actually yields the best fit. This function was 
also shown to yield an excellent fit to the Terminal Veloc- 
ity Curve of the Milky Way Galaxy [1] , and was found to 
work better than the standard function to fit the disper- 
sion profile of the globular cluster system of NGC 4636 
in the transition regime between the Newtonian and the 
deep-MOND regime js^. Even more interesting is to 
check that the stellar M/L ratios implied by this simple 
/i-function are realistic. On average, we have: 



(M,/L)e,^i 
(M,/L) 



= 0.70 ±0.15. 



standard 



(9) 



i.e. their value is ~30% lower than when using the stan- 
dard MOND interpolating function. In TE\ the standard 
MOND M/L ratios are slightly higher than the predic- 
tions of stellar population synthesis models (by ~ 0.10 
dex, using his fits of the stellar M/L vs colour plots for 
B-V>0.55, i.e. above the break), for a given observed 
colour. Recently, [l^l compared the ranges of M/L ra- 
tios allowed by different methods to determine them, and 
suggested that the normalisation of the stellar M/L vs 
colour relation should be lowered by 0.05 to 0.10 dex. 
Note that they consider the predictions of the models by 
[ilj , while [3| uses the most recent models by (i^ . How- 
ever, in the range of colours considered by ^5], 0.4 < B-V 
< 0.8, the two models give predictions that differ by at 
most 0.03 dex. The stellar M/L ratios that we find from 
the simple /i-function, ~ 30% lower than the standard 
MOND ones, fall perfectly in this range allowed by the 
various constraints considered by [40| . 

Let us end this section by noting that the "traditional" 
view on this observational baryon-gravity relation would 
be that it is a purely empirical relation that will be under- 
stood within the framework of the concordance ACDM 



cosmological model. Even though this is a very challeng- 
ing task, the meaning of the a-family of Eq. ^ combined 
with Eq. (1) would then correspond (in spherical symme- 
try) to the following empirical relation: 



9Nb 

9dm 



agNb + gpM 
ao 



(10) 



where gNb and goM are the contribution to the rotation 
curve gravity from the baryons and DM respectively. If 
neither gravity nor inertia are to be modified, the present 
success of the simple relation a = 1 to fit galaxy rota- 
tion curves still captures the essence of the fine-tuning 
problem for dark matter. The mass-to-light ratios of Ta- 
ble I being compatible with current population synthesis 
models, Eq. pO)) tells us that the ratio of baryonic grav- 
ity to DM gravity should be equal to the total gravity 
in units of oq. For example, if one plots the contribu- 
tion of baryons and dark matter to the rotation curve at 
a certain radius as a function of the baryonic contribu- 
tion to gravity at that radius (Fig. 4), one finds a cross 
reminiscent of Il5l. 



III. THE ^i-FUNCTION AND THE FREE 
FUNCTION OF TeVeS 

A. TeVeS in galaxies 

Here, we present the theoretical constraints upon the 
/i-function in the framework of Bekenstein's relativistic 
multifield theory of gravity 25j. Before doing that, we 
need to briefly look back on the basics of the theory, and 
recap the relevant equations (in Sects. Ill A and III B). 

If one wants to modify Newtonian gravity on galaxy 
scales in order to get MOND, one can modify Poisson's 
equation in order to recover Milgrom's law (second equal- 
ity in Eq. 1) in highly symmetric systems |2J], namely: 



V.[/i(|V$|/ao)V$] -47rGArp, 



(11) 



where Gat is the usual Newtonian gravitational constant, 
<I> is the MOND gravitational potential and p is the den- 
sity of the (baryonic) matter distribution. In spherical 
symmetry, one then immediately finds back Milgrom's 
law (Eq. 1) from Gauss's theorem. Out of spherical sym- 
metry, one finds that 



/i(|V<I>|/ao)V$ = V$Ar + V X h, 



(12) 



where h is a regular vector field determined (up to a gra- 
dient) by the condition that the curl of the total force V'l' 
must vanish. Inside spiral galaxies, a numerical MOND 
potential solver Q, |4^ must thus be used to obtain 
the exact MONDian force field corresponding to Eq. (|lip . 
However, it has been shown that Milgrom's law is still a 
good enough approximation to fit galactic rotation curves 
within that framework [H, |4^ . 

On the other hand, the extension of General Relativity 
that has been proposed by Bekenstein for the MOND 
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FIG. 2: Rotation curve decompositions with a — 1 [see Eq.(|4} and Eq.(l2}] and ao — 1.35 x 10~* cm s~^. The distances were 
left as free parameters, within 30% of the adopted distances in the original papers. The solid line is the best fit to the data, 
the dotted/dashed/long-dashed lines represent the disk/gas/bulge Newtonian contributions. Note that the Newtonian curves 
deviate from the MOND ones even at small radii, which is characteristic of the simple function a — 1. 
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FIG. 3: Rotation curve decompositions with the standard MOND /^-function [see Eq.(|3])], and ao = 1.21 x 10~* cm s~^. The 
distances were left as free parameters, within 30% of the adopted distances in the original papers. The solid line is the best fit to 
the data, the dotted/dashed/long-dashed lines represent the disk/gas/bulge Newtonian contributions. Note that the distances 
were already allowed to vary in fits with the standard function in pjTj ] (and that the fits presented on this figure are thus the 
same) but not in i9l. 
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FIG. 4; This figure illustrates that Eq. dTOjl captures the 
essence of the fine-tuning problem for CDM in galaxies. For 
a = 1, we plot the square root of the fraction of total gravity 
accounted by the baryons (dashed line) and the DM (full line) 
as a function of the logarithm of baryonic gravity. The cross 
is reminiscent of the one obtained in p^ . fig. 4]. 

paradigm is a bi-metric multi-field theory ^25*1 , namely a 
tensor-vector-scalar (TeVeS) theory leading to Milgrom's 
law in spherical symmetry as we shall see below, but 
different from Eq. (jlip in more general geometries [ssj . 
Interestingly, this theory can be reformulated in various 
fashions, such as e.g. a pure tensor-vector theory in the 
matter frame [l^]. We however concentrate here on the 
original formulation of the theory. 

The tensor of the theory is the Einstein metric gafs out 
of which is built the usual Einstein-Hilbert action, while 
the MONDian dynamics comes from a dynamical scalar 
field (p whose lagrangian density is aquadratic. All the 
matter fields are coupled to a physical metric involving 
the tensor and the scalar field of the theory [2^ eq.(21)], 
and which must be disformally related to the Einstein 
metric in order to consistently enhance the deflection of 
photons. This needs the introduction of a dynamical nor- 
malized vector field Ua ■ The total action is thus the sum 
of the Einstein-Hilbert action for the Einstein metric, 
the action of the scalar field, the action of the vector 
field, and the matter action involving the physical met- 
ric. Einstein-like equations are then obtained for each 
of these fields by varying the action w.r.t. each of them 

The action of the scalar field 4> can be written 

where G is the bare gravitational constant appearing in 
the Einstein-Hilbert action, g is the determinant of the 



Einstein metric, and y = c~^kPh°'f^da(f>dp(j) (the metric 
h'^^ being some specific combination of the Einstein met- 
ric and the vector field [1^ ) . Here fc is a small dimension- 
less parameter (of the order of one percent) determining 
the strength of the scalar coupling to matter, I is a length 
scale linked to the acceleration scale ao, and c is the speed 
of light. The function /(y) is a free function driving the 
dynamics of the scalar field, that can be identified (to a 
factor) with its lagrangian density [ill, and that yields 
the MOND behaviour in galaxies when judiciously cho- 
sen. Note that for quasi-static systems (for which time 
derivatives are zero) such as galaxies, the scalar field is 
akin to k-essence scalar fields, that were notably intro- 
duced as possible dark energy fluids that could also drive 
inflation [il, |4^, [s^l . This name comes from the fact that 
their dynamics is dominated by their kinetic term — /(?/), 
contrary to other quintessence models in which the scalar 
field potential plays the crucial role. 

In the Einstein frame, matter fields couple to the 
three gravitational fields (tensor, vector and scalar field) 
through an effective "disformal" metric that differs from 
the Einstein metric, and which therefore, although dy- 
namical, assumes a prior geometric form. In that frame 
the above k-essence field may thus be considered to me- 
diate a "fifth force" gs = — V(/), in addition to the usual 
gravity connected to the gravitational tensor. On the 
other hand, shifting our viewpoint from the Einstein 
frame to the matter one, matter fields appear to couple 
to some metric that assumes no prior form, but which is 
determined by a non-diagonal system of coupled differ- 
ential equations. In that case the scalar field does not 
trivially manifest as a fifth force. Those two points of 
view are however equivalent, and we shall now work in 
the Einstein frame. 

If one couples such a k-essence field to the matter sec- 
tor in the way it is done in TeVeS, one finds in the non- 
relativistic limit, neglecting pressure, and in a static con- 
figuration, the following modified Poisson's equation (25l . 
eq.(42)]: 

V. [Wcj^f (c-4fcZ2|V0|2)] = kGp, (14) 

where p is the density of the matter distribution [25[. 
The analogy with Eq. (|lip is striking. For a spherical 
distribution of mass M one thus gets 

Wcf>f{c-^kfm') = ^ (15) 

outside the matter. It must be stressed that many pa- 
pers achieve to find modified theories of gravity in which 
a test particle has an asymptotic circular velocity of the 
form V — where 7 is a pure number of the theory, 
and thus seem to recover fiat rotation curves of spiral 
galaxies [Sll, [sl, . However one must not only be con- 
cerned with the existence of the plateau velocity, but 
also with the whole phenomenology of rotation curves 
(see Sect. II and e.g. 113, [E^), including the TuUy-Fisher 
relation [H, [s^ . This relation teaches us that 7 should 
behave as M^^'^ (at roughly a ten percent level due to 
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the possible difference between tlie asymptotic velocity 
and the velocity at the last observed radius). In these 
theories it is not correct to set by hands such a value 
for 7. On the contrary, k-essence theories achieve this 
in a very elegant manner. It is indeed immediate to see 
that if /'(y) behaves as then the "fifth force" domi- 
nates the Newton- Einstein gravity in the ultra-weak field 
limit, and both the asymptotic plateau velocity and its 
amplitude are recovered. K-essence theories would not 
have been useful if, say, the TuUy-Fisher relation were 
(X M with p ^ 4. In other words, k-essence theo- 
ries are relevant toy-models in this context because the 
effective (MOND) gravitational field in the outskirts of 
galaxies is the square root of the Newtonian one. 



B. TeVeS and Milgrom's law 

To the leading order, the physical metric near a quasi- 
static galaxy in TeVeS is given by the same metric as 
in General Relativity, with the standard Schwarzschild 
potential replaced by the total potential 



(16) 



where (/)„ is the Newtonian potential obtained from using 
the bare gravitational constant G, and S 1 is a dimen- 
sionle ss p arameter which slowly evolves with the cosmic 
time [23, eq.(58)]. Note that this parameter can always 
be chosen to be precisely S = 1 at the present time. This 
means that, locally and at the classical level, the scalar 
field (f) plays exactly the role of the dark matter poten- 
tial. To the leading order, we have the following field 
equations in a quasi-static system: 



(17) 



i.e. the Poisson equation for the Newtonian potential, 
and 



V.[M,(|V0|/ao)V</)] =4^Gp, 



(18) 



i.e. a MOND-like equation (see Eq. [TT|) for the scalar 
field, obtained from Eq. by defining (noticing that 
2/ > in a static configuration) 



(19) 



Note that the acceleration scale oq does not need to be 
precisely constant through cosmic time. However, since 
the variation is extremely slow, it can be considered as 
effectively constant during the entire time that galaxies 
have existed (see Sect. Ill C). From Eqs. ([TH), (HT]), and 
(fT8| we are now able to look for a Milgrom's law for 
the total gravitational potential Note that we restrict 
ourselves to spherical symmetry so that we look for an 
equation of the form 



/i(|V$|/ao)V$ 



GNM{r)r 



(20) 



together with the usual asymptotic MONDian conditions 
/i(a;) ~ a; if a; <C 1 and /^(x) ~ 1 if a; ^ 1. The above 
equation notably implies that Gat is the gravitational 
constant measured in Newtonian systems (and hence on 
Earth) , which is not required to be equal to the bare one. 
For commodity we define 



ao 



ao ' 



|V^ 

ao 



9_ 

a-o 



(21) 



the MONDian gravitational force per unit mass and its 
modulus in units of ao, and 



ao 



ao ' 



ao 



£1 

ao 



(22) 



the scalar field force and its modulus in units of ao. Given 
that ao can slowly (even though not much) vary with 
cosmic time , the very definition of s slowly varies too. 
In spherical symmetry, we thus have 



1 



fJ-s{s)s = — (x - s) 



|V0» 
ao 



(23) 



Let us emphasize that the fis function can be quite arbi- 
trarily chosen in a k-essence theory (but see Sect. IIV!^ . 
so that the /z-function that will be derived from it may 
not possess the right MONDian behaviour. We shall 
therefore adopt a phenomenological point of view her- 
after, by imposing the standard asymptotic conditions 
for /i, and look for the corresponding properties of /is. 

Let us define a parameter vq that connects the renor- 
malized gravitational constant measured on Earth to its 
bare value: Gjv = i^qG. In spherically symmetric TeVeS, 
Milgrom's law can then be written as 



fi{x)x = -;=r{x — s), 
meaning that the /i-function is defined by 



1 + ^Ms 



(24) 



(25) 



Note that in [2^, sect. IV B] Milgrom's function was de- 
fined by 



jl{x)x 



ao 



(26) 



neglecting the renormalization of G as a first order ap- 
proximation. Since it is always possible to immediately 
renormalize G by a factor S such that the total potential 
in Eq. (fT6)) is the sum of the Newtonian potential and the 
scalar field, it is tempting to define jj. = "E^jl. However, 
as we shall show hereafter, this renormalization is only 
valid for some special forms of the TeVeS free function. 
More generally the exact //-function of Eq. (p5)) is related 
to the one of Eq. (pS)) by /i = vof^i and the fact that vo 
may differ from S has wide implications. 

In the rest of this section, we exhaustively examine how 
the free function of TeVeS, hereafter /is(s) as defined in 
Eq. ([in]), connects with the /x-function of MOND, and 
we insist on the importance of this famous value of vo 
inherited from the shape of /i^. 
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C. The deep-MOND regime 

In the deep-MOND regime, we must have ^i{x) ~ x for 
X ^ 1, which imphes, from Eq. that s ^ 1 (s ~ a;) 
and from Eq. (gSl) that 



s 



(27) 



in this regime. In order to get this behaviour in Eq. ([19] 
with the standard choice f'{y) ~ \/vT^ for small y [1^ 
we must define for the length scale: 



3k 



(28) 



Note that since I is a parameter of the theory (fixed by 
definition), the above equation shows that v^a^ must be 
a constant. As we will see below however, the renormal- 
ization of G, encoded in the parameter vq depends on the 
cosmic time, since z/q is related to S. This implies that 
the acceleration scale ag does vary with cosmic time (ag is 
higher at higher redshift). Let us underline that what we 
call ao is the quantity determining the MONDian asymp- 
totic circular velocity for a spherical system, since in the 
decp-MOND regime test-particles feel a force given by 
^fGMU^/r by virtue of Eqs. (gSl), dM]) and (H?]). This 
force is interpreted by an observer on Earth as a force 
given by ^/GJjMmo / r , where Gat = vqC is our local con- 
stant of gravity. However Bekenstein [25;] estimated that 
the time variation of S is extremely mild in the matter 
dominated period, so that as far as fits of galaxy rotation 
curves are concerned, this variation of ag is expected to 
have a negligible influence, even though ag is expected to 
be somewhat different at the epoch of nucleosynthesis. 

Note also that Eq. (^5)1 agrees with [H, eq.(62)] only 
in the very special case where = 2. This special case 
is also very important in the study of the Newtonian 
regime, as we shall see hereafter. 



D. The Newtonian regime: unbounded /^s-function 

The case I'o = S is a physically interesting one, be- 
cause, a priori, it is always possible to immediately renor- 
malize G in Eq. p?)) . such that E!(/)„ = $jv. Then, also 
replacing G by Gn = SG in the scalar field modified 
Poisson equation (Eq. fTS]) means that the phenomcnolog- 
ically relevant interpolating function for the scalar field 
becomes S/Zs (during a period for which the parameter S 
is constant, which is effectively the case during the entire 
time that galaxies have existed), and that the particular 
case 1^0 = 2 thus corresponds to no a posteriori renor- 
malization of the gravitational constant. Eq. [23]) and 
Eq. ([M]) can then be rewritten: 



s = ^{x)x ~ ^/is(s)s 



(29) 



The question is then whether it is actually possible 
to have = S. Let us first show that we indeed have 



fQ = S if and only if the /i^-function is unbounded, i.e. 
fis oo. Indeed, if vq = S, then, from the first equality 
in Eq. ((29|) , we can write /i as 1 — s/x where we must have 
s <ti X for the Newtonian regime a; 3> 1, which proves 
from the second equality in Eq. (|29p s — x/ {l + Ens) that 
fis must diverge with s. Respectively, if /is(s) diverges 
for some value s = sg, it is straightforward from Eq. (j25|) 
that /i — > I'o/S and therefore that i/g = S by definition 

of fl. 

Then let us emphasize that Eq. ((29l) actually corre- 
sponds to the implicit assumptions in 5(|, where it was 
shown that, without any a posteriori renormalization of 
the gravitational constant, it was impossible to find a 
TeVeS free function /is(s) corresponding to the standard 
/i-function of MOND (Eq. [3|) . This particular result can 
actually be generalized through the following theorem: 
When i/g = S, a MOND /i-function can be obtained from 
a TeVeS free function only if n < 1 in its asymptotic 
expansion /i(a:) ~ 1 — A/x^^ (where A is some positive 
constant). 

The proof is quite straightforward: 

Since = 'E., we have from Eq. (j24p (or from the first 
equality in Eq. 



a;[l - fi{x)] 



(30) 



Let us assume that n > 1, then for a; — s- oo we would 
have s ~ A/x"^^ 0. But since s ^ 1 is already the 
deep-MOND regime (see Sect. Ill C), it cannot be also 
the Newtonian one □ 

What it means is that, as can be easily checked [s^ . 
the /is function would be multivalued. This is however 
not acceptable since /is must define unambiguously the 
dynamics of the scalar field. No unbounded /is-function 
can thus reproduce the standard /i-function (Eq. [3]) , for 
the reason that the latter admits the following asymp- 
totic expansion 



/i(x) - 1 - 



1 

2^' 



(31) 



and by virtue of the previous theorem (here n — 2). 

On the other hand, if n < 1, this is not a problem since 
s would diverge jointly with x, while if n = 1, s would 
asymptote to a finite positive value, and the function 
s{x) from Eq. pO]) could still be strictly increasing, and 
thus be invertible. E.g., the simple /i-function of MOND, 
(Eq. [2|), has n — 1 for its asymptotic expansion. Its 
corresponding function s(x) from Eq. (j30p is 



s{x) 



1 



leading to the unbounded /is function 
S/is(s) 



1 - s 



(32) 



(33) 



This corresponds to the n — case of the free function 
proposed by [H, eq.(13)] for the local value S = 1. The /i- 
functions tested against observations in Sect. II also have 
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an unbounded /i^ counterpart in TeVeS. If we want their 
corresponding k-essence free functions to be expressed as 
pure functions of the constants in the action of TeVeS 
(see Sect. Ill A), independently of cosmic variations of 
S, we can write the a-family of Eq. (|4]) as: 



(34) 



Another interesting aspect of these unbounded fig func- 
tions (or of the special case fo = S) is that they lead to 
an analogy with electromagnetism that might be help- 
ful to grasp the feeling of the to-be-understood physical 
meaning of the scalar field of TeVeS. Indeed, with the 
definition of Eq. (|25p . one can then write the Newtonian 
gravity generated by a baryonic point mass M as: 

GpfM fi 

— o— ^ gN ^ [gN + gs) = x .g^- (35) 

In general, the conservative electric field E is given by 

D P 

E = d + p, with d EE — , p = , (36) 

Co Co 

where eo is the absolute permissivity of the vaccum, D is 
the electric displacement vector, and P is the polarization 
vector. If we consider a dielectric medium of relative 
permissivity /i surrounding a free charge ~Q, we have 



— = d^fix{d + p) 



1 



A* 



X p. 



(37) 



The analogy with TeVeS is striking, the free charge —Q 
playing the role of a baryonic point mass M (both creat- 
ing a l/r^ force in the absence of the "dielectric"). The 
scalar field " fifth" force gs plays the role of the polariza- 
tion field p and the total gravitional acceleration g plays 
the role of the conservative electric field E. The electro- 
magnetic dielectric permissivity fi is often constant for a 
uniform medium, depending on the microscopic structure 
of the medium, but can vary spatially in the presence of a 
non-uniform field. Interestingly the polarization field in 
a polarizablc dielectric medium has also an microscopic 
couterpart: the randomly oriented dipoles in the medium 
are lined up to form a macroscopic dipole (of amplitude 
P per unit volume) when the electric field E is applied. 
Could this hint towards a quantum origin for a TeVeS- 
like field theory, with the vacuum playing the role of the 
dielectric medium, and the scalar field force arising from 
the gravitational effect of the baryonic matter? This is 
all too speculative for the time being, and this analogy 
should of course break down at the cosmological level 
in TeVeS because of the role played by the vector field. 
This is pretty much beyond the scope of this paper, but 
we hope the analogy may help the reader to grasp the 
feeling of one possible concrete physical meaning of the 
scalar field. 

What is more, it is actually striking that, for the best- 
fit a = 1 function of Sect. II and at S = 1, the ^-factor 



of Eq. (1), is nothing else than the scalar field strength 
itself in units of ao, i.e. 



(38) 



from Eq. ^ and Eq. (|32p . This may acquire a great 
theoritical significance in any future TeVeS-like theory 
where the /i-function could be viewed as a dielectric-like 
factor. The fact that what could be interpreted as a "di- 
gravitational" relative permissivity of space-time could 
simply be given by the modulus of the scalar field gra- 
dient is very appealing and may pave the way towards 
a well-motivated relativistic theory of MOND. We stress 
that this relation only holds locally in TeVeS because 
of the time-evolution of S (although this parameter is 
effectively constant during the entire matter dominated 
period), but that this factor comes from the vector field 
in TeVeS, and could thus perhaps be avoided in some- 
what different relativistic theories of MOND [23l[57l. 



E. The Newtonian regime: bounded ^s-function 

We are now going to show that the above theorem 
constraining the MOND /i-function in the context of un- 
bounded TeVeS /ij-functions breaks down when i^q differs 
from S (something which was not realized in [56j). 

If 1^0 7^ S, then, using the condition fj.{x) 1 for 
X ^ oo (i.e. the MOND condition to recover Newtonian 
dynamics in the high acceleration regime), we find from 
Eq. (gll) and Eq. that s > 1 and that the Hs function 
is bounded, i.e. asymptotes to a finite value which reads 



1 



Mo = 



1^0 -:^ 



(39) 



Note that fig must be strictly positive (see Sect. IIV A[) 
so that I'o > S. Respectively, it is straightforward from 
Eq. (|25p that if /is(s) — Mo for s — > oo, then we have that 
fi foMo/(l + S^o) = 1, meaning that i/q = S -I- l/no- 

The above theorem for unbounded functions then 
breaks down because, in order to get a /i-function with 
an asymptotic expansion /i(x) ~ 1 — A/x", we need a 
scalar field strength such that 



n+l 



1 + ^/io 1 + ^^J.o 



(40) 



meaning that, when x is large, s does not tend to zero, 
and fJ.sis) is not necessarily multivalued when n > 1. 

It is actually possible to find such bounded fig func- 
tions that give rise to the standard /i-function of Eq. ([3]). 
There is in fact a one-parameter family of such functions 
(labelled by the value fis{oo) — /io, or equivalcntly by 
vq). When inserting the standard /i-function of Eq. ^ 
in Eq. (|24p . one must make sure that the function 



1 - —n{x) 
fo 



(41) 
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is strictly increasing, and thus inve rtible. This means 
that me must have I'o/S > 4/3-^/2/3. 

However, such bounded /Ltg-functions lead to a non- 
trivial renormalization of the gravitational constant, and 
their relation with the /^-function of MOND is much less 
elegant than in the case of unbounded functions, leading 
to a breakdown of the analogy with clcctromagnetism. 
Indeed, in the case of unbounded functions, the relevant 
function for the scalar field could be considered to be S/i^ , 
and one needed S/x^ = s for small s to recover fj,{x) = x 
for small x. With the definition of I in [2^, eq.(62)], this is 
exactly the behaviour of /is for Bekenstein's toy model at 
small s. However, since Bekenstein's /Xg is a bounded iunc- 
tion, its corresponding MOND /^-function is proportional 
to j/qx/S instead of x in the weak regime, because of the 
renormalization of G. This means that, in order to re- 
cover the TuUy-Fisher law in weak gravity, I must be rede- 
fined as in Eq. ([28| . But the /i- function of Eq. ([5]) is then 
only a rough approximation of Bekenstein's toy model, 
even in the intermediate regime, because of the renormal- 
ization of G that makes the /^-function asymptote more 
quickly to /i = 1 than Eq. On the other hand, note 
that there exists a one-parameter family of bounded /is- 
functions precisely yielding Eq. ([5]), but different from 
Bekenstein's free function, while for a given S the un- 
bounded function yielding Eq. ([5]) is simply S/i^ s. 



F. Non-Spherical systems 

The relations between fi and /is we derived above are 
of course only valid in spherical symmetry. Out of spher- 
ical symmetry, Milgrom's law in Eq. (1) is not exact for 
any MOND-like gravity theory: the Newtonian force, the 
classical MOND force, and the scalar field force are no 
longer parallel. The curl field obtained when solving for 
the full classical MOND potential $ in Eq. ((TT]) wiU be 
different from the one obtained when solving the equa- 
tion for the scalar field (f> (Eg. [T5|) (see e.g. _35:]). Note 
that, for an unbounded /is-function {i^q = S), the latter 
can be rigorously remoulded into the following expression 



S = Sn + gs 



V • [Efisgs - Sn] = 0. (42) 



Here the total gravity g consists of a baryonic Newto- 
nian gravity g^r (calculated from the Newtonian Pois- 
son equation with the Newtonian gravitational constant 
Gn — SG) plus gs as an extra or boosted gravity (the 
space part of the four vector gradient of the scalar field cj) 
of TeVeS , simply playing the role of dark matter gravity) . 
The amount of boost is specified by S/is, which takes the 
following dual meaning: (1) the ratio of baryonic gravity 
w.r.t. extra-gravity, or (2) the ratio of effective inertial 
mass moff = (not to be confused with the actual in- 
ertial mass, unmodified in a modified gravity theory) to 
the effective loss of inertial mass m(l — /i), with fi defined 
as in Eq. (^5)1 . 

Here, for the first time in the literature, we compare 
the different computed rotation curves for the same disk 
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FIG. 5: Bottom: exact Newtonian (dot-dashed line), classical 
MOND (dashed line), and TeVeS (full line) equatorial-plane 
rotation curves for a disk galaxy model with density distribu- 
tion given in equation (|43|l . The dotted line is the correspond- 
ing rotation curve approximated using Eq.(l). The top panel 
is a zoom-in of the upper left corner of the bottom panel. 



galaxy using the prescriptions of Eq. (1) (Milgrom's law), 
Eq. dlT]) (classical MOND), and Eq. ^ (or equivalently 
Eq. U^l i.e. TeVeS). We consider an exponential disk 
model typical of a high surface brightness galaxy, char- 
acterized by a baryonic density distribution 



Pb{R,z) = 



47ri?o ^0 



exp 



'#)(cosh^)-^ (43) 



with total mass Mi, — 4.2 x lO^°A/0, and characteristic 
scale lengths i?o = 3 kpc and zq = 0.1 kpc, where R and 
z are the cylindrical coordinates. 
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We can compute the exact force field, by solving 
Eq. pT|) using the numerical MOND potential solver de- 
veloped by [ij]. This solver was originally tested in [3| 
for the standard /i- function of Eq. ([3]), but we verified 
it works as well in solving Eq. (fTTj) for other /i-functions 
such as Eq. ^ and Eq. (O, or in solving the scalar field 
Poisson equation for a /x^-function such as Eq. (|33p . 

For this model we compute exactly four different 
equatorial-plane rotation curves, which are plotted in 
Fig. [51 the Newtonian rotation curve (dot-dashed) , the 
classical MOND rotation curve for the simple /i-function 
of Eq. (121) (dashed) , and the TeVeS rotation curve for the 
/is-function of Eq. ([55)) (full line; we assumed for this plot 
ao = 1.2 X 10-8 cm s^^). The MOND/TeVeS rotation 
curves deviate significantly from the Newtonian one even 
at small radii, which is characteristic of the simple func- 
tion (see Sect. II). Fig. [Slalso plots as a dotted fine the 
approximated rotation curve calculated from Eq. (1) with 
the simple /x-function. From the bottom panel of Fig.[5|it 
is apparent how for an exponential disk model the exact 
TeVeS rotation curve, the exact classical MOND rota- 
tion curve, and the corresponding approximated rotation 
curve are practically indistinguishable from an observa- 
tional point of view. Only a zoom- in at intermediate radii 
(Fig. [3 top) allows to appreciate the difference between 
the three modified-gravity curves: the TeVeS curve is 
systematically slightly higher (by less than 0.5%) than 
the classical MOND one, while the approximated curve 
using Milgrom's law slightly overestimates the velocity 
(by only a few tenth of percent) for R < 2Ro, and slightly 
underestimates it (by less than 1%) for R > 2i?o- Note 
that if we consider i? = 25 kpc as the typical last ob- 
served radius in this high surface brightness galaxy, the 
difference between the circular velocity at that point and 
the asymptotic one (~ 5-10 km/s) is in any case higher 
than the difference between the velocities predicted at 
that point from Milgrom's law and from TeVeS due to 
the non-sphericity of the system (~ 1-2 km/s). 



IV. COSMOLOGY 
A. The Cauchy problem 

We have seen that the asymptotic condition ^{x) ~ x 
when X <^ I implies that /is(s) must behave linearly with 
s when s is small (Eq. l27p . and notably that we must have 
/is(0) = 0. This s = signals the transition between local 
physics and cosmology in TeVeS [25j . 

Let us however emphasize that a k-essence theory such 
as TeVeS can exhibit superluminal propagations when- 
ever f"{y) > [11]. Although it does not threaten 

causality [53|, one has to check that the Cauchy prob- 
lem is still well-posed for the field equations. Within the 
choice of signature (— , 4-, + , +) of local Lorentzian met- 
rics, it has been shown [HI, HI, HOl that it requires the 
otherwise free function / to satisfy the following proper- 



ties, \/y: 

f'iv) > (44) 
f'{y) + 2yf"{y) > 0. (45) 

In terms of the scalar field strength s = {c? / la^) y / k 
(when J/ > 0, for a quasi-static galaxy), and of the auxil- 
iary function /Zs(s) = 47r/'(j/)/A; , these conditions read, 
Vs > 0: 

^ls{s) > (46) 
Ai.(s)-f s/4(s) >0. (47) 

where the prime denotes the derivative with respect to s. 
Note that, throughout Sect. Ill, although we have shown 
that it was sometimes impossible to recover a given /i- 
function from any single- valued unbounded /ig-function, 
we implicitly assumed that for a given /i^ it was always 
possible to recover a spherically symmetric formulation 
"a la Milgrom" of local physics of TeVeS, i.e. an equation 
like Eq. (j20p . This means that when inserting a given 
/is in Eq. ([23[) . the inversion of x{s) was assumed to be 
always possible (i.e. x{s) is strictly increasing). Making 
use of S > 0, this is actually guaranteed by the second 
inequality above, meaning that sfis{s) must be a strictly 
increasing function of s. 

Now, concerning the transition between local physics 
and cosmology in TeVeS, the asymptotic behaviour of 
the /i-function can therefore not be implemented by a 
consistent k-essence theory, since /is (0) = would violate 
both inequalities at s = 0. As explained in [s^, Issl ]. 
all present day k-essence-like theories of MOND such as 
TeVeS predict the existence of a singular surface around 
each galaxies on which the scalar degree of freedom does 
not propagate, and can therefore not provide a consistent 
picture of collapsed matter embedded into a cosmological 
background. 

A simple solution consists in assuming a modified 
asymptotic behaviour of the /i-function, namely of the 
form 

fi{x) r^eo + x (48) 

if a: <C 1. In that case there is a return to a Newtonian 
behaviour at a very low acceleration scale a; <C eoj E^nd 
rotation curves of galaxies are only approximatively flat 
until the galactocentric radius 

1 IGnM 

Rr J^I . (49) 

eo V ao 

One must thus have eg <C 1. A similar phenomenology 
also arises in a slightly different context in [6l[. Present 
data on galaxy rotation curves ||52] suggest that l/£o 
must be at least of the order of 10. Note that Eq is a 
quantity which slowly evolves with cosmic time (such as 
and ao). In order to have such a behaviour for /i, 
one must have in the deep-MOND regime /is ^ s/vq + e 
for small s, and we then have — vqE (where e is a 
constant). 
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Note that, in the negative range of values of y corre- 
sponding to the time-like sector of the expanding Uni- 
verse, can still be defined as a function of y as in 
Eq. fTO]) . It can also be viewed as a function of s if we 
extend the definition of s to 

s^sygn{y)^M. (50) 

Then the Cauchy conditions read, Vs < 0: 

Ms(s) > (51) 
^ls{s) - s^i'M >Q. (52) 

B. The mirror-image function 

To be able to describe both galaxies and cosmology, 
the free function of TeVeS must thus be able to deal also 
with non-stationary systems such as the Universe itself. 
We describe hereafter how this can be achieved, following 

Bekenstein's original proposal [25| was to construct the 
scalar field action of Eq. (fT^ whith a lagrangian density 
given as a one-to-one function of /is viewed as an auxiliary 
non-dynamical scalar field, i.e. f{y) = J-{^s)- 

Such one-to-one construction had the drawback that 
the lagrangian density of the scalar field necessarily had 
unphysical "gaps", i.e. that a sector was reserved for 
space-like systems (e.g. from dwarf galaxies to the solar 
system in the range Q < Hs < for y > 0) and a dis- 
connected sector was reserved for time- like systems (e.g., 
the expanding Universe in the range /is > 2/io for y < 0). 
While viable mathematically, such a disconnected Uni- 
verse would not permit galaxies to collapse continuously 
out of the Hubble expansion. 

In an effort to re-connect galaxies with the expanding 
Universe in TeVeS, it was proposed in [56| to construct 
the lagrangian density as a one-to-one continuous func- 
tion of y, rather than a function of /is, thus not consid- 
ering jjLa as an auxiliary non-dynamical scalar field any- 
more. This allows for a smooth transition from the edge 
of galaxies where to the Hubble expansion. It was 
also suggested in [5g| to extrapolate the lagrangian den- 
sity for galaxies into the cosmological regime by simple 
mirror- imaging, to minimize any fine-tuning in TeVeS, 
i.e. 

/is(s) = /is(-s) (53) 

with s defined as in Eq. ([50]) . meaning that the time-like 
sector is a simple mirror-image of the space-like sector. 

In other words, it means that, in order to produce 
a consistent cosmological picture, we suggest to extend 
the a-family of Eq. to the time-like sector by using 
the following k-essence function f{y) in the action of the 
scalar field fEq. [T3)) : 

f{y)^f{0) + l''^dy (54) 



with 

fis = e+ ^ , < £ < 0.1, < a < 1, (55) 

Ao - ay/\y\ 

where Aq = V3fc/47r. 

It would thus be of great interest to compute the CMB 
anisotropics and matter power spectrum [H, [l^ using 
this action for the scalar field allowing for a smooth tran- 
sition from the edge of galaxies to the Hubble expansion. 
Note also that the integration constant /(O) can play the 
role of the cosmological constant [6^ (/(O) — corre- 
sponding to no cosmological constant), but that even a 
no-cosmological constant model could drive late-time ac- 
celeration |66|, [B^I , which is not surprising since k-essence 
scalar fields were first introduced to address the dark en- 
ergy problem. This TeVeS cosmological model is thus 
worth exploring since it could have the ability to resolve 
many problems at the same time. However, the a-family 
of /i-functions does not work very well on galaxy cluster 
scales if gjv is the gravity generated by baryons alone 
in Eq. This problem of the MOND-phenomenology 

has been known for a while [i^ , and it was then suggested 
[6^ that ordinary neutrinos of 2 eV could be included 
in the "known matter-total gravity" relation. Including 
gAT as the gravity of both baryons and neutrinos, it was 
found in [TOj] that a reasonable fit to the lensing of the 
now famous bullet cluster 0] is possible if neutrinos have 
a 2 eV mass (a mass also invoked in current MOND fits of 
the CMB anisotropies [H,!?!!), while the unusually high 
encounter speed of the bullet cluster [t^ could be repro- 
duced since the modified gravity potentials are deeper 
than the usual CDM ones. This cosmological framework 
was called the /xHDM paradigm in [t^] by contrast with 
the usual ACDM one. These neutrinos would not change 
the gravity on galaxy scales due to their low densities 
[lOI, hence would not affect the goodness of the rota- 
tion curves fits presented in Sect. II. However, for ellip- 
tical galaxies at the very center of galaxy clusters where 
the neutrino density is the highest, the velocity disper- 
sion profiles could be affected by these neutrinos, thus 
explaining the poor purely baryonic fits for NGC 1399 
when using Eq. ([2]) [7^ (even though the external field 
effect in this galaxy, due to the gravity of the cluster it- 
self, might render the fits difficult to achieve in practice 
even when including neutrinos). 



V. CONCLUSIONS 

Throughout this paper we have studied from a theoret- 
ical point of view the observationally motivated baryon- 
gravity relation in galaxies [14, 15, 16]. This relation is 
completely encapsulated within theu-function (see Eq. 1) 
of Modified Newtonian Dynamics [l3| , each distinct func- 
tion defining a distinct baryon-gravity relation. This rela- 
tion might be empirical or fundamental, but the present- 
day success of the relation with the sole constraints of 
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the MOND paradigm, i.e. /i(a;) = 1 at x ^ 1 and 
^{x) = X at X <C 1 with x ~ g/cHo, points at first sight 
to a modification of gravity in the ultra- weak field limit. 
Even though the most important aspect of MOND re- 
mains this asymptotic behaviour, the study of the transi- 
tion between the Newtonian regime and the deep-MOND 
regime, characterized by the shape of the /i-function, may 
also be important in the quest for a possible fundamental 
theory underpinning the MOND prescription. We have 
thus investigated the observational constraints upon the 
/i-function in that transition zone, as well as the theoret- 
ical constraints upon it in the framework of the relativis- 
tic formulation of MOND dubbed TeVeS Note that 
the latter are only valid if MOND really has its basis in a 
multifield theory like TeVeS, which is of course far from 
guaranteed. 
We have shown that: 

i) Observationally (from the quality of galaxy rotation 
curves and the current accuracy of population synthe- 
sis models), one cannot distinguish between the "sim- 
ple" (Eq. [ID and the "standard" (Eq. O form of the 
/z-function, but the transition from Newton to MOND 
cannot be much more gradual than that implied by the 
simple function (Sect. II). 

ii) The simple form of /i does require generally lower 
mass-to-light ratios for the stellar components of spiral 
galaxies than does the standard form. This could, com- 
bined with accurate population synthesis models, even- 
tually distinguish between them (Sect. II). 

iii) A particular class of TeVeS free functions (the un- 
bounded fis functions in Eq. [TH] and Eq. correspond 
to a trivial renormalization of the gravitational constant 
Gjv = SG where S is a parameter of the theory, and are 
easily linked to the /^-function of MOND. However, this 
/i-function must have n < 1 in its asymptotic expansion 
/i(a;) ~ 1 — A/x" (where A is some positive constant), 
meaning that the standard form of fi (Eq.[3]) cannot be re- 
covered from a single- valued unbounded fig (Sect. Ill D). 

iv) Renormalizing the gravitational constant as Gat = 
vqC with j/Q > S for bounded /is functions allows to 
recover /i-functions with a steeper asymptotic expansion. 
The standard /i of Eq. ^ can be recovered if is such 
that Eq. dU) is invertible (Sect. Ill E). 

v) Observationally (from galaxy rotation curves) one 
cannot distinguish between Milgrom's law [l^, classical 



MOND ^ and TeVeS [H. Using the code of [3 we 
have shown that the difference between the circular veloc- 
ity at the last observed radius in a typical spiral galaxy 
and the asymptotic circular velocity is higher than the 
difference between the velocities predicted at that point 
from classical MOND and from TeVeS due to the non- 
sphericity of the system (Sect. Ill F). 

vi) Some modification of the /i-function is necessary to 
make cosmology consistent with quasi-static mass dis- 
tributions in TeVeS, which should lead to a return to 
1/r^ attraction at very low accelerations in galaxies 
(Sect. IV A). 

vii) By extrapolating the scalar field lagrangian den- 
sity for galaxies into the cosmological regime by simple 
mirror- imaging in TeVeS, one minimizes the fine-tuning 
of the theory while avoiding that the sector reserved for 
space-like systems such as galaxies be disconnected from 
the sector reserved for time-like systems such as the Uni- 
verse itself (Sect. IV B). 

Let us finally note that the simple form of /i (Eq. ^ , 
yielding excellent fits to galaxy rotation curves and not 
needing a non-trivial renormalization of the gravitational 
constant in TeVeS, has to be distorded in the limit of 
high accelerations to be in accordance with observed 
planetary motions in the inner solar system [36]. This 
could be achieved for an unbounded /ts by e.g. somehow 
making the parameter a dynamical in Eq. (|34p . going 
smoothly from 1 in galaxies towards 20 at the orbit of 
Mercury. More prosaically, it can be achieved by mak- 
ing the function bounded as in the n = 1 case of [56l . 
eq.(13)]. 
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